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Abstract 

A system with two correlated Gaussian white noises is analysed. This system can describe 
both stochastic localization and long tails in the stationary distribution. Correlations be- 
tween the noises can lead to a nonmonotonic behaviour of the variance as function of the 
intensity of one of the noises and to a stochastic resonance. A method for improving the 
transmission of external periodic signal by tuning parameters of the system discussed in 
this paper is proposed. 

Key words: Correlated Gaussian noises, stochastic resonance 
PACS: 05.40.Ca 



1 Introduction 



Stochastic models with multiplicative, or parametrical, noise find numerous appli- 
cations in a variety of branches of science and technology. Unfortunately, mod- 
els for which analytical results are known are very scarce and any such a model 
deserves a thorough discussion. Recently, Denisov and Horsthemke [1] have dis- 
cussed a model given by the equation 



x = —ax + \x\ a r](t) , (1) 



where ^ a ^ 1, r](t) is a Gaussian noise, possibly coloured, and have found 
that it can describe anomalous diffusion and stochastic localization. Denisov and 
Horsthemke have also discussed several physical systems in which models of the 
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type (1) can be useful; see references provided in their paper. Later Vitrenko [2] 
has generalized (1) to include two noise terms: 



x = — 



(a + m(t))x + \x\ a r] 2 (t) , 



(2) 



where rji^ are certain coloured and correlated Gaussian noises. This system has 
a very nice feature: for < a < 1, it interpolates between a linear transmitter 
with multiplicative and additive noises (a = 0) and a linear system with a purely 
multiplicative noise (a = 1). The two limiting cases, a = and a = 1, are 
very well known in the literature (see e.g. Ref.[3] and references quoted therein). 
Vitrenko has formally linearized the system (2) by means of a substitution that has 
been already used in [1]: 



and solved the resulting equation for the trajectories. Converting back to the orig- 
inal variable proves to be rather tricky and that author has managed to do so only 
if the noises rj 1:2 are correlated in a very specific (not to say peculiar) manner. It is 
now widely recognized that correlations between various noises can lead to many 
interesting effects. It is, however, possible that phenomena reported by Vitrenko re- 
sult principally from the very specific form of correlations assumed by this author 
and are not generic to the system (2). We find it interesting to see how the system 
behaves for the intermediate values of a when the correlation requirements are less 
restrictive than those discussed by Vitrenko. 

Coloured noises introduce more complexity. However, if a dynamical effect is 
present in the white noise case, it also appears, perhaps in a distorted form, in the 
coloured case [4a-b]. To simplify the discussion, we will assume that the noises 
are white. Finally, note that the expression a + r)i(t) in Eq. (2) can be interpreted as 
a biased noise. The noise that multiplies \x\ a in Eq. (2) is not biased. To "symme- 
terize" the system, we include a bias in £ 2 m our analysis. It is also convenient to 
have explicit expressions for noise amplitudes, or coupling constants between the 
noises and the dynamical variable. We thus recast the equation (2) in the form 



where a > 0, ^ a ^ 1, £i )2 are mutually correlated Gaussian white noises: 



x = — 



(a + p^(t))x+\x\ a (b + qUt)), 



(4) 



&(*))= 0, (fc(f)fc(f )> = *(*" , < = 1,2, 



(5a) 
(5b) 
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and c G [—1, 1]. If not otherwise specified, we interpret the noises in the sense of 
Ito. For the sake of terminology, we will call the noise "multiplicative" and 
£ 2 (t) "additive", even though this terminology is accurate only if a = (for a > 
both noises couple parametrically). Note that if a particle hits x — 0, it stays there 
forever if a > 0. 

There is, in fact, one more reason for including b ^ in our discussion. Much as 
the substitution (3) linearizes the system (4), another substitution, namely 



z = £L (6) 

x 



converts it to a noisy logistic equation 



z = (1 - a)(a + P Ut))z - (1 - a)(6 + g&(*))* 2 • (7) 

We have discussed this last system in [6] and found that 6^0 together with cor- 
relations between the noises can lead to a nonmonotonic behaviour of the variance 
(z 2 ) — (z) 2 as a function of the intensity of the "additive" noise, q, and to a stochas- 
tic resonance [7] if the system is additionally stimulated by an external periodic 
signal. It would be naive to expect that these phenomena occur in the system (4) in 
exactly the same manner as they do in (7). A nonlinear change of variables, espe- 
cially in case of stochastic equations, can significantly alter the behaviour. We will 
see, however, that there are striking similarities between the systems (7) and (4). 

This paper is organized as follows: We construct the Fokker-Planck equation for 
the system (4) in Section 2 and in Section 3 we present its stationary solutions. 
Then in Section 4 we discuss the constructive effects of the correlations between 
the noises; in particular, in Section 4.2 we give numerical evidence for the presence 
of the stochastic resonance. Conclusions are given in Section 5. 



2 The Fokker-Planck equation 



The problem of constructing a Fokker-Planck equation corresponding to a process 
driven by two correlated Gaussian white noises has been first discussed in Ref. [8], 
where the two noises have been decomposed into two independent processes. The 
same result has been later re-derived in [9], where the authors have attempted to 
avoid an explicit decomposition of the noises but eventually resorted to a disguised 
form of the decomposition. The Fokker-Planck equation for correlated white noises 
has been also discussed in Refs. [10] and [11] and in several other papers; see, for 
example, Ref. [6] for a particularly simple re-derivation. 
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A general Langevin equation 



x = h(x) + gi(x)^(t) + g 2 (x)Ut) , (8) 

where x(t) is a one-dimensional process and £ lj2 are as in Eqns. (5), leads to the 
following Fokker-Planck equation in the Ito interpretation: 



d -^T 1 = -^(*)P(*> + \^ B ^) P ^ t) , (9a) 

where 

B(x) = [ 9l {x)} 2 + 2c gi (x)g 2 (x) + [g 2 {x)] 2 . (9b) 
In case of Eq. (4) we obtain 



dP(x,t) d, 1 d 



- — (ax-b\x\ a )P(x,t) + - — (p 2 x 2 - 2cpqx\x\ a + q 2 \x\ 2a ) P(x,t) . 
at ox I ox v J 

(10) 



In the following we interpret \x\ a as \x\ a = (x 2 ) a ^ 2 , where the square of x must be 
calculated prior to taking the fractional power. It is also apparent that the probability 
P(x, t) does not depend on the absolute signs of the amplitudes p, q, but only 
on their relative sign. We assume that sgn(pg) = +1. This comes at no loss of 
generality as the equation (10) is invariant under a simultaneous change of signs of 
pq and c. 

Finding stationary distributions corresponding to Eq. (10) is the main goal of this 
paper. This, in principle, could be handled by standard methods [12a-b], but it 
would be very difficult due to the absolute value and the fractional powers. It is 
apparent that since the right-hand-side of the corresponding stationary equation 
vanishes identically if x — 0, the term S(x) should always be included in any 
stationary distribution. We now use the substitution (3). After some algebra we 
eventually obtain 



ay-b + ^-{p 2 y 2 - 2cpqy + q 2 ) 



dt dy 



P(y,t) 



+ -(1 - a) 2 — (p 2 y 2 - 2cpqy + q 2 )P(y, t) . (11) 



The last term in the square brackets in Eq. (11) corresponds to the Ito interpre- 
tation [13]. This term is missing if the noises are interpreted according to Stra- 
tonovich. 
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The stationary distribution solves an equation that is fairly easy to integrate: 



(a+(l- ^aj p 2 ^jy-b- cpq + ° q 



2y 



p, 



St 



1 /i _ \ /_ 2 2 o | 2 \ ^~st 



+ -(1 - a)( P Y - 2cpay + q 2 )-^ = . (12) 
2 ay 



3 Stationary distributions 

Before proceeding to the general case, let us discuss the case where there is only 
one noise present. 

3.1 No "multiplicative" noise 



If there is no "multiplicative" noise, p — 0, and no bias in the "additive" noise, 
6 = 0, our problem reduces to that discussed in Ref. [1]. Eq. (12) takes the form 



Kfh+^-^iH- <i3) 

This equation corresponds to the following Langevin equation 



y = -(ay+ ^ +VI- aqUt) , (14) 
which, in turn, corresponds to an overdamped motion in a potential 
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Vcs(y) = -ay 2 + -aq 2 In \y\ . (15) 

The effective potential (15) has an infinite noise-created well at y = which 
traps Brownian particles; this well is missing if the noises are interpreted accord- 
ing to Stratonovich. Curiously, in another context we have observed a similar phe- 
nomenon, where noise interpreted according to Ito created an insurmountable bar- 
rier restricting particles to one half of the real axis [14]. A similar barrier is ob- 
served in the noisy logistic system (7), cf. Ref. [6]. Loosely speaking, the change 
of variables (6) converts an infinite barrier into an infinite well. 
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The equation (13) is solved by 



N ( ay 2 \ 

PAy) = ]yf^ eXP \^a)i^ (16) 

where N is a normalization constant, or transforming back to the original variable 



I I a ex P /-, \ 

\x\ V (1 — a)q z 



P tt (x) = Tns exp\- ^ '— ■ (17) 



The distribution (17) is normalizable for a > and ^ a < 1. For a = it 
reduces to a standard Gaussian distribution, and for < a < 1 it mildly diverges at 
x — 0, where some particles are trapped, or stochastically localized. If a increases 
towards unity, the divergence becomes more pronounced, or more particles become 
localized. At the same time, though, tails of the distribution get heavier which is 
characteristic for anomalous diffusion. If we interpret the system (4) with b = as 
"interpolating" between linear systems with an additive and multiplicative noises, 
we can see that if a = 0, the stationary distribution is non-singular. The "less 
additive" the system becomes as a increases, the more pronounced the singularity 
is and the tails of the distribution get flatter. Eventually, for a linear and purely 
multiplicative system, the distribution reduces to 8(x) as all particles collapse to 
the origin of the force. 

The presence of a bias, b ^ 0, introduces some asymmetry in the exponential term, 
but the overall behaviour remains much the same: 



N (2bx\x\~ a -a(x 2 f " ' 

i 7a ex P 7\ \ n 

\x\ \ (1 — a)q z 



P*(x) = ^ex P [ n ^ . (18) 



If the noises are interpreted according to Stratonovich, we obtain 



and there is no stochastic localization. 



3.2 No "additive" noise 



If there is no "additive" noise, q — 0, and no bias, 6 = 0, the only normalizable 
stationary solution is P st (x) = S(x), corresponding to all particles eventually col- 
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lapsing to their common resting point. If b ^ 0, there is no stationary solution as 
some particles go to the resting point, but some can escape to infinity. 



3.3 The general case 

If the noises are not maximally correlated, |c| 7^ 1, the general solution reads 



where N is again a normalization constant. Despite its complicated form, principal 
properties of the distribution (20) are easy to find. Because the inverse tangens 
function, arctan(-), is limited, the exponential term is also limited and convergence 
properties of (20) depend solely on its denominator. One can easily see that this 
distribution is normalizable for all ^ a < 1. For < a < 1 the stochastic 
localization, in the Ito interpretation, occurs. The distribution has rather heavy tails. 
It has a convergent first moment if a > \ap 2 . The second moment is convergent if 
a stronger condition, a > |(1 + a)p 2 , is satisfied. 

If either b 7^ 0, or c 7^ 0, or both, the distribution (20) is not symmetric. Apart from 
the narrow singularity around x — 0, it has another peak centred on the minimum 
of q 2 — 2cpqx\x\~ a + p 2 (x 2 ) 1 ~ a . Its location depends on the sign of c: if c > 0, the 
peak is located to the right of x = 0, and if c < 0, it is located to the left. If |c| < 1, 
the height of this peak can be very large. Thus, in the stationary state, a majority of 
Brownian particles is stochastically localized either around the central singularity, 
or in the additional peak created by the correlations. The tails do not contribute 
much to the overall density. However, if a approaches unity, the tails, and the tail 
with the same sign as the location of the peak in particular, become rather heavy 
and outliers, or particles far removed from both the peak and the singularity, can 
easily be found. The asymmetry between the tails is introduced by the exponential 
term: The distribution is further asymmetrically broadened by the exponential term 
in the numerator of (20). This broadening can be removed if 



It is important to understand the origin of this phenomenon. The asymmetric broad- 
ening results from the bias — the force acting in one direction is, on the average, 
larger than the force acting in the opposite one. In the system (4) the parameter 
b 7^ acts as one source of the bias; it has been introduced for this specific pur- 
pose. It is also known that correlations between two noises can effectively introduce 
another bias, see e.g. Refs. [3,8,10]. If the condition (21) is met, the two sources 




(20) 



bp — caq = . 



(21) 
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of bias nullify each other. To see this, let us represent the two correlated Gaussian 
white noises £ lj2 as linear combinations of two independent GWNs ^ lj2 : 



&(t)= ^i(t), 

( 2 (t) = c^(i) + ^^W. 



(22a) 
(22b) 



With the condition (21) satisfied, the Langevin equation (4) now takes the form 



x 



{a + pMt)) (x-^x^ +VT^q\x\ a tp 2 (t). (23) 



The system now behaves as if it were driven by two uncorrelated white noises, one 
of which is unbiased. As a result, the bias-induced asymmetric broadening of the 
stationary distribution disappears. 



3.4 Maximally correlated noises 



The distribution (20) does not have a universal limit |c| — > 1. Instead, if c = ±1, 
we need to solve Eq. (12) directly and then convert back to the original variable. 
We obtain a candidate solution 



Atrial (^) 



eX P((T^ 



2(b P Taq) 



)p 2 (qTpx\x\ a ) 



\X\ 



q =F px\x\ 



2 ( 1+ (T^) ' 



(24) 



where the =p sign is the opposite of the sign of the correlation coefficient, c = ±1. 
However, the right-hand- side of (24) is not normalizable. If q =f px\x\~ a = 0, the 
exponential in (24) hits its essential singularity. This singularity is eliminated if 
a special case of the condition (21), namely 



bp =F o,q = , 



(25) 



holds. In this case, either p = q = and the system becomes fully deterministic, or 
the stationary Fokker- Planck equation (12) factorises: 



2py 



[(2a + (2 - a)p 2 )y =f apq]P st (y) + (1 - a)py(py T q) 



dPM 
dy 



0. 
(26) 



8 



A singular distribution 5(py =F <?) solves Eq. (26). The regular part of this equation, 
the one in the square brackets, again leads to a not normalizable solution. We, 
therefore, conclude that if the noises are maximally correlated and the condition 
(25) holds, the stationary distribution reads 

P st (x) = 7 S(x) + (1 - 7 )5 (x T (q/p) 1/(1 - a) ) , (27) 

where 7 is the fraction of the initial population that collapses to zero; observe that 
with our sign convention adopted, q/p > 0. If the noises are maximally correlated 
but the condition (25) is not satisfied, there is no stationary distribution. There is 
a striking similarity between the system (4) and (7), where a similar situation oc- 
curs [6]: If a condition analogous to (25) is satisfied and the noises are maximally 
correlated, the noisy logistic system has a 5-like stationary distribution. If the noises 
are maximally correlated but the counterpart of the condition (25) does not hold, 
a stationary distribution does not form in the noisy logistic system, either. 



4 Constructive effects of correlations 



As we have seen, a delicate interplay between the correlations and the bias can 
significantly alter the shape of the stationary distribution. We may expect that this 
can lead to various unexpected properties of the system (4). 



4. 1 Nonmonotonic behaviour of the variance 

Recall that depending on the parameters, the stationary distribution of the system 
discussed in this paper can be nearly limited to very narrow peaks; with a different 
set of parameters, these peaks can be asymmetrically broadened. The second central 
moment of a probability distribution, (x 2 ) — (x), if convergent, is perhaps one of its 
simplest and most easily comprehended characteristics. It is interesting to see how 
the second moment of the distribution (20) behaves as a function of the "additive" 
noise strength. Because of the complicated analytical structure of this distribution, 
we have not been able to evaluate the integrals ff^ x P st (x) dx, ff^ x 2 P st (x) dx 
analytically. We have done so numerically instead. Example results are presented 
in Fig. 1; parameters chosen correspond to a convergent second moment. As we 
can see, a minimum of the variance as a function of the "additive" noise strength 
is clearly visible. This minimum is fairly deep if the correlations are large and 
becomes very shallow as the correlations decrease. Note that if b < 0, the minimum 
appears for negative values of the correlation coefficient (not plotted). 
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0.250 0.375 0.500 0.625 0.750 0.875 1.000 



q 

Fig. 1. Nonmonotonic behaviour of the variance of the distribution (20) as a function of the 
"additive" noise strength, q, for various correlations between the noises, c. Other parameters 
are a = 1/8, a = 1.25, b = 1.0, p = 1.0. 

4.2 Stochastic resonance 

Now suppose that the system discussed in this paper is additionally stimulated by 
an external, periodic signal. The Langevin equation takes the form 

x = -(a+p^(t))x+ \x\ a {b + qi 2 (t) + A cos(Qt + </))) (28) 

where the noises are as in (5). Because we do not know exact solutions of a time- 
dependent Fokker-Planck equation corresponding to Eq. (28), we have solved the 
equation (28) numerically with the Euler-Maryuama algorithm and a timestep equal 
2~ 12 . To generate the correlated noises £ lj2 , we have first generated two indepen- 
dent Gaussian white noises ^ lj2 ; we have used the Marsaglia algorithm [15a-c] for 
that purpose and the famous Mersenne Twister [16] has been used as the underlying 
uniform generator. Then the correlated noises are created as linear combinations of 
the two uncorrelated ones, see Eq. (22) above. Example trajectories of the system 
(28) and associated power spectra, averaged over 128 realizations of the noise and 
on the initial phase of the signal, <fi, are presented in Fig. 2. The shape of the trajec- 
tories and the power spectra strongly depend on the parameters of the system, and 
on the correlation coefficient, c, and the strength of the "additive" noise, q, in par- 
ticular. In general, the higher the correlations, the more ordered the trajectories are. 
It is worth noting that higher harmonics of the driving frequency can be visible in 
the power spectra, indicating a nonlinear nature of the coupling between the signal 
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Fig. 2. Panel (a): a fragment of a typical realization of the process (28) with c = 1. Panel 
(b): the corresponding power spectrum averaged over 128 realizations. Panels (c), (d): same 
as (a), (b) above, but with c = 0.5. Other parameters, common for all panels: a = 1/8, 
a = 1.25, b = 1.0, p = 1.0, q = 0.8, A = 1, and = 2tt. 

and the dynamical variable. 

To quantify these observations, we will use the Signal-To-Noise Ratio (SNR) as 
a measure of the stochastic resonance: 



SNR = 10 log 



5. 
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signal 



Snoise(/ = Q/271) ' 



(29) 



where Signal is the height of the peak in the power spectrum at the driving frequency 
and 5' no i se (/) is the frequency-dependent noise-induced background. Several other 
measures of the stochastic resonance have been proposed [17a-c], but we choose 
the SNR as the simplest, oldest and most commonly used one. Selected results, 
averaged on both realizations of the noises and the initial phase, are presented in 
Figs. 3 and 4. For high values of the correlation coefficient, a clear maximum in the 
SNR is visible. This shows that there is an optimal level of the "additive" noise that 
maximizes the ratio of power transmitted through coherent oscillations induced by 
the driving signal to that transmitted by the irregular ones, or that there is a stochas- 
tic resonance in the system (28). For correlations only slightly larger than zero, the 
resonance is very small and it disappears for c ^ 0. Note that this happens if the 
asymmetry parameter, b, is greater than zero. For b < the stochastic resonance 
occurs for negative correlations and reaches its largest magnitude at c = — 1 . In the 
symmetric case, 6 = 0, there is no stochastic resonance. Again, these features of 
the stochastic resonance resemble very much those of the noisy logistic system (7) 
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Fig. 3. Signal-to-noise ratio for the system (28). Parameters are a = 1/8, a = 1.25, 
6 =1.0, p = 1.0, A = 1, and = 2tt. Curves presented correspond, back to front, to 
the following values of the correlation coefficient: c = 1.000, 0.875, 0.750, 0.625, 0.500, 
0.375, 0.250, 0.125, and 0.000. 




Fig. 4. Same as Fig. 3 but with a = 7/8. 

discussed in [6]. 

The resonance becomes sharper if a approaches unity (Fig. 4). At the same time, 
values of SNR away from the resonance are much smaller than those in the small 
a case. 
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Fig. 5. Signal-to-noise ratio for the system (28) as a function of the bias parameter, b, and 
the correlation coefficient, c. Other parameters are a = 1/8, a = 1.25, p = 1.0, q = 0.8, 
A = 1, and n = 2ir. 

4.3 Response to a change of deterministic parameters 



We have shown in the two previous Subsections that the system discussed here can 
be optimised by choosing an appropriate level of the "additive" noise. In practice, 
however, controlling the amplitude of the noise or correlations between the two 
sources of the noise can be very difficult. Tuning the deterministic part of the system 
may be much easier to achieve, and as our discussion of the asymmetric broadening 
of the distribution in Subsection 3.3 has shown, by changing the bias parameter, 
b, we can optimise the system even if the noise amplitudes and the correlation 
coefficient are not known. 

To test for that, we have again numerically simulated the externally stimulated sys- 
tem (28) by the same means that have been used in Subsection 4.2 above. This time 
amplitudes of the two noises have been kept constant and the bias parameter has 
been varied. Selected results are presented in Fig. 5. As we can see, changing the 
bias does optimise the system. Clear maxima in the signal-to-noise ratio are visible. 
These maxima are most pronounced if correlations are large, |c| < 1, but they are 
present also for |c| ~ 0, even though the overall shape of the curves is much flatter. 
For the uncorrelated case, c = 0, the weak maximum coincides with b = which is 
to be expected due to symmetry of the system. To put it in a slightly different way, 
we can see that the uncorrelated system transmits an external signal badly. Any cor- 
relations between the noises potentially improve the transmission properties. The 
system can be optimised to reach its full potential by appropriately adjusting its 
deterministic parameters. 
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5 Conclusions 



In this paper we have discussed a nonlinear system with two correlated sources of 
Gaussian white noises. A closely related system has been discussed previously by 
Vitrenko in Ref. [2]. We have been mainly interested in what happens when the 
restrictions on correlations between the noises imposed by that author are lifted 
and, additionally, when the "additive" noise becomes biased. We have shown that 
this system can display both stochastic localization and heavy tails in its stationary 
distribution which is characteristic for anomalous diffusion. This agrees with pre- 
viously published results [1,2]. It is worth noting, though, that authors of that Ref- 
erences obtained their results under the assumption that the noises were coloured; 
we have shown that the same happens for white noises as well. 

Next, we have shown that correlations present in the system discussed here can lead 
to interesting constructive effects of the noise: to a nonmonotonic behaviour of the 
variance of the stationary distribution and to a stochastic resonance. Finally, we 
have shown that the system can be optimised to an external periodic signal not only 
by varying amplitudes of the noises, but also by tuning the deterministic parameters 
of the system when the noise amplitudes and the correlation coefficient between the 
noises remain, in principle, unknown. 

Surprisingly, the system (4) discussed here is related to the noisy logistic system 
(7) that we have discussed previously [6]. As we have shown, many, but not all, 
properties of these two systems are strikingly similar. 
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